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Abstract 

We introduce and investigate by numerical simulations a number of 
models of emotional agents at the square lattice. Our models describe 
the most general features of emotions such as the spontaneous emotional 
arousal, emotional relaxation, and transfers of emotions between different 
agents. Group emotions in the considered models are periodically fluctu- 
ating between two opposite valency levels and as result the mean value of 
such group emotions is zero. The oscillations amplitude depends strongly 
on probability p s of the individual spontaneous arousal. For small values 
of relaxation times r we observed a stochastic resonance, i.e. the signal to 
noise ratio SNR is maximal for a non-zero p a parameter. The amplitude 
increases with the probability p of local affective interactions while the 
mean oscillations period increases with the relaxation time r and is only 
weakly dependent on other system parameters. Presence of emotional 
antenna can enhance positive or negative emotions and for the optimal 
transition probability the antenna can change agents emotions at longer 
distances. The stochastic resonance was also observed for the influence of 
emotions on task execution efficiency. 

1 Introduction 

Recently physicist are interested in modeling various social phenomena, e.g. 
opinion evolution, culture migration or language dynamics [TJ [H [31 IH [Z] ■ 
One of such new domains of investigations is emotion research O [9] [10] . Emo- 
tions or affective states can be caused by external or internal processes but dif- 
ferences between these two emotion sources are sometimes hard to distinguish 
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[TTj . An external emotion-arousing event comes usually as a result of inter- 
actions with environment. The case when emotions are changing or emerging 
without any influence of external factors corresponds to spontaneous emotional 
arousals. The important feature of an affective state is its short life time [TTj . 
Phenomena of emergence and decline of an emotion are usually fast and unex- 
pected processes which are frequently difficult to control or to predict. In other 
words: emotions evolve on the spur of the moment. In the situation when nei- 
ther external nor internal emotion generating processes take place, the emotion 
relaxes to a non-affective (neutral) state. In general, one can say that emotions 
are positive or negative, i.e. they possess a positive or a negative valence. What 
is important, emotions are always directed towards somebody or something [T2"j . 
It means that the emotional valence has an influence on actions of an individual 
and can also affect other people. As a result, in a moment of affective interac- 
tion we can define an emotional emitter and an emotional receiver. It is obvious 
that they can change their roles in time. 

Social interactions take place in social space [T3]. As we observe in every- 
day life, people are usually much more influenced by emotions of their rela- 
tives, friends, acquaintances or partners with common goals than by emotions 
of strangers. Social relations or common goals frequently increase emotional 
influence [5]. It follows that social distances should be introduced to properly 
describe affective phenomena. 

In our paper we introduce and investigate by numerical simulations a number 
of models of emotional agents. The models take into account such features as 
spontaneous emotional arousal, emotional relaxation and affective interactions 
between agents. We study also the influence of emotions on tasks execution 
efficiency. 

2 Model description 
2.1 Model 1 

We consider the behaviour of agents at the two-dimensional square lattice 
with the unit lattice constant and the size X ■ Y. The agent density is p = 1 
thus the total number of agents is N = X ■ Y. The agents do not move in space 
but the agent i can change his emotional state (t) = +1, 0, —I in the course of 
time. For numerical simplicity we shall assume that the time variable is discrete 
t = 1,2,3.... The state e, (t) = will be called emotionally neutral, while the 
state ej (t) = I (—1) is emotionally positive (negative). An emotional state of 
an agent evolves as a result of three processes: due to inter-agents emotional 
interactions, due to a process of spontaneous emotional arousal, and due to an 
emotional relaxation. 

The square lattice represents the social space and the social distance between 
agents i, j will be assumed as the smallest number of edges at a path between i 
and j, i.e. rij = \xi — Xj | + \yi — y 3 ■ |, where Xj, yi describe Cartesian coordinates 
of the agent i. 
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Table 1: Scheme of interactions between emotional emitter and receiver in Model 
1 



Initial emotional states of agents ei(0) are randomly selected with the uni- 
form probability distribution 1/3. 

During the simulation, changes of emotional states are observed. As men- 
tioned before, we distinguish the following processes leading to changes of emo- 
tional states ej(f). 

(i) Affective interactions between agents are collective effects. We randomly 
select an agent i that will be considered as an emotional emitter. Since 
people usually interact only with somebody whom they are closely related 
to, emotions of the emitter can influence with the probability p all agents j 
(emotional receivers) that are placed at a distance ry < e from the emitter 
i where e is the range of emotional interactions. The emotional emitter i 
can change an affective state of emotional receiver j following the Table 
ffl 

(u) The spontaneous emotional arousal means in our model a transition from 
the current emotional state to any state (including the current one) with 
the same probability p s /3. It follows the process corresponds to a kind of 
internal emotional noise. 

(iii) The emotional relaxation process is described as follows. If the emotion 
e, of the agent i and emotions of all agents in his closest neighbourhood 
(r = 1) do not change over time window r, then the emotional state of the 
agent relaxes to the non-emotional value, e, = 0. Similarly to the process 
of spontaneous emotional arousal, every agent is considered separately in 
this scheme. 

The whole simulation consists of T time steps. In one step we consider first 
the emotional relaxation effect, then the spontaneous emotional arousal (for all 
agents) and at the end affective interactions between agents. In every time step 
we randomly select ./V emotional emitters influencing their neighbourhood. 
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2.2 Model 2 



Our Model 1 suffers from many simplifications. To improve it, we increased 
heterogeneity of agents interactions in the Model 2. Individuals more likely 
interact with close people. For this reason for any agent i we divide agents 
into three groups: friends, acquaintances and strangers. Agents i, j are friends 
when Tij — 1, i.e. they are closest neighbours. The distance r.y = 2 defines 
acquaintances, while r,j > 2 means strangers. 

Likewise in the Model 1, agents can change their emotional state as a re- 
sult of relaxation, spontaneous emotional arousal and affective interactions. The 
difference is the way agents interact with each other. Friends interact with prob- 
ability p, and this is the most likely interaction. Acquaintances have a smaller 
influence on each other and we assumed the probability of their interactions is 
p/2. There are no interactions between strangers. 

3 Results of numerical simulations 

3.1 Group emotion 

All presented simulation results were received for the same size of the system 
X = Y = 40. We are interested in the group behaviour during the time of sim- 
ulation. For this purpose, we examine a group emotion as an average emotional 
state of the group at time t: 



Although the updating rules for positive and negative emotions in our model 
are symmetric (no kind of valence is favored) , one could suppose that collective 
interactions between agents would lead to the spontaneous symmetry breaking 
as for example in the two-dimensional Ising model [14]. We found that the 
spontaneous ordering was absent and the group emotion averaged over a long- 
time of simulation was always close to zero even when the value of the noise 
parameter p s was very small. Instead of the spontaneous order there are large 
amplitude oscillations of the group emotion (e) (t) in such a limit. We stress that 
the oscillations exist although the system dynamics does not contain any incrtial 
part. We suppose that the oscillations are a result of individual relaxation 
processes in the considered models. Fig. 1-6. present the time dependence of 
(e) (t) for selected values of p, p s and r and one can see the influence of these 
parameters on a typical frequency of group emotion oscillations and on their 
amplitudes. 

The dependence of the oscillations amplitude on the spontaneous arousal 
probability p s is non-monotonic, see Fig|I][2j and the amplitude is maximal 
for a characteristic value p s « 0.1 (when one fixes other system parameters as 
p = 1 and r = 2). The increase of the transfer probability p leads to a monotonic 
increase of oscillation amplitude, see Figj3]|4] Values of p and p s possess only a 
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Figure 1: Behaviour of group emotion (e) (t) for Model 1 with r = 2, p = 1 and 
for different values of p s = 0.01, 0.1, 0.6, 1 in A-D, respectively. 
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Figure 2: Behaviour of group emotion (e) (t) for Model 2 with r = 2, p = 1 and 
for different values of p s = 0.01, 0.1, 0.6, 1 in A-D, respectively. 
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Figure 3: Behaviour of group emotion (e) (t) for Model 1 with r = 2, p s = 0.1 
and for different values of p = 0.1, 0.25, 0.5, 1 in A-D, respectively. 
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Figure 4: Behaviour of group emotion (e) (t) for Model 2 with r = 2, p s = 0.1 
and for different values of p = 0.1, 0.25, 0.5, 1 in A-D, respectively. 



8 





Figure 5: Behaviour of group emotion (e) (t) for Model 1 with p s — 0.1, p — 1 
and for different values of relaxation time t — 2,3, 5, 10 in A-D, respectively. 



weak influence on the period of observed oscillations. The typical period of these 
oscillations increases with the relaxation time r, see Fig[5] and [U Comparing 
the Model 1 and the Model 2 for the same parameter values we observe larger 
amplitude oscillations in the first one. 

For a quantitative analysis of these oscillations we have performed a Fourier 
transform 

X k = ^(e)(t)exp{~kt). (2) 
t=o 

Some representative plots of the amplitude of the Fourier components 
are presented at Fig. [7] where one can observe a characteristic peak in the 
frequency domain. The peak decays for large values of the noise parameter 
p s . Let us introduce a parameter reminding the signal-to- noise ratio (SNR) 
used for observations of stochastic phenomena [TS] that describes the relative 
strength of this peak and is calculated as: 

SNR= f^nax\ (3) 

where A max is the height of the peak in the Fourier transform while A av is the 
mean value of this transform (averaged over all frequencies). Plots of SNR as 
a function of p s show a clear maximum for some intermediate value of p s , see 
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Figure 6: Behaviour of group emotion (e) (t) for Model 2 with p s = 0.1, p = 1 
and for different values of relaxation time r = 2, 3, 5, 10 in A-D, respectively. 
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Figure 7: Amplitude of Fourier transform of group emotion (e) (t) for Model 1 
with T = 1024, p = 1 and r = 3 for different values of p s = 0.01, 0.15, 0.6, 1 in 
A-D, respectively. 



Fig. EJ Such a behaviour reminds the well known phenomenon of stochastic 
resonance [THl HE1 HZ] that occurs in a large class of dynamical systems where 
a system nonlinearity interferes with external noise. In our system a typical 
resonance behaviour with a single maximum exists when the relaxation time 
r is not to large and for large r the SNR can possess several maxima, see 
plots in Fig. 03 and [5p. The inverse of the characteristic frequency of the 
maximum peak in the Fourier spectrum increases with the relaxation time r 
until a saturation effect occurs for large r (see Fig(Hl). This frequency increases 
slightly with the parameter p s when the parameter r is small (see Fig |10|) . The 
influence of the interaction probability p on this frequency is very weak. 

For small values of p s we observed forming a spatial group of agents who 
possessed the same emotional state and as result there were periods of time 
when the whole group was emotionally polarized. This effect was enhanced 
when the transfer probability p was close to 1 (each emotional emitter interacts 
with emotional receivers in a deterministic way) and large collective oscillations 
in mean group emotions were observed in such a case. 
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B) Model 1 for t=5 and t=1 
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C) Model 2 for t=2 and x=3 
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D) Model 2 for t=5 and x=1 
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Figure 8: Signal-to- noise ratio (SNR) as a function of p s with p — 1 for different 
values of r = 2, 3, 5, 10 for Model 1 (A-B) and Model 2 (C-D). 
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B) Model 2 




Figure 9: Characteristic time (inverse of characteristic frequency) as a function 
of t for p — 1 and for different values of p s for Model 1 and Model 2 in A-B, 
respectively. 
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A) Model 1 
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Figure 10: Characteristic frequency as a function of p s for p = 1 and for different 
values of t = 2, 3, 5, 10 for Model 1 and Model 2 in A-B, respectively. 
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Table 2: Probabilities of interactions for a positive antenna in Model 1. 
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Table 3: Probabilities of interactions for a positive antenna in Model 2. 



3.2 Emotional antenna 

Up to now we have assumed that the ability to transfer own negative/positive 
emotions is the same for every member of the considered social group. We have 
assumed also that the ability to percept and to share other people emotions is 
the same for every agent. 

In reality both these features are related to individual affective character- 
istics that vary from human to human and are also dependent on specific in- 
teractions links [5]. To model this effect we introduce a simple heterogeneity 
to the considered agent model assuming that some agents are more likely to 
transfer their own positive or negative emotions and the same agents are also 
more likely to share positive or negative emotions of other agents. We call such 
type of agents emotional antennas and we distinguish between positive and neg- 
ative emotional antennas. By definition a positive (negative) emotional antenna 
favours to transfer its positive (negative) emotion to surrounding agents and it 
is also very likely to be influenced by positive (negative) emotions of neighbour- 
ing agents. Although the antenna possesses a different capability of interactions 
with other agents (depending on its emotional state) it can be in the positive, 
negative or neutral state and it is governed by the same rules of the emotional 
relaxation and spontaneous arousal as other agents. 

Interactions probabilities between a positive antenna and its neighbourhood 
are presented in Tables [2] and |3l For the negative antenna the bottom row 
corresponding to e a = — 1 should be exchanged with the top row corresponding 
to e a = 1. Of course one can also imagine other types of emotional antennas. 

During computer simulations we located a positive antenna in the central 
position of the lattice and we observed how its presence influences emotions of 
neighbouring and more distant agents. At Fig. [TT] we observe that the mean 
emotion of agents at the distance larger than 2 is nearly zero. The exception is 
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Figure 11: Average emotional state of antenna's neigbourhood as a function 
of distance between antenna and neighbours r (Model 1 and Model 2 in A-B, 
respectively) for r = 2, p s = 0.05 and different values of p. 
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A) Model 1 B) Model 2 
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Figure 12: Average emotional state of antenna's neigbourhood as a function 
p for different values of distance between antenna and neighbours r for r = 
p s = 0.05 (Model 1 and Model 2 in A-B, respectively). 
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Figure 13: Average emotional state of antenna's closest neigbours (r = 1) as a 
function of p s for different values of p for r = 2 (Model 1 and Model 2 in A-B, 
respectively). 
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Figure 14: Number NC of successful agents over a time T — 1024 with upper 
potential barrier U u — 30 and lower potential barrier Ud = — 50 as a function of 
p and p s for Model 1 and r = 2. 

the case p s» 0.2 where the influence of the antenna is felt at larger distances, 
see Fig. [121 The result can be understood as follows. When p = 1 the positive 
antenna broadcasts its negative emotions with a probability 0.5 thus its positive 
role is not very significant as compared to the case of lower values of p. However 
when p << 1 then the positive antenna does not distribute negative emotions 
but it also does not receive positive emotions from its neighbourhood which is 
mostly in the neutral state = 0. It follows there is an optimal value of the 
transfer probability p when the role of antenna is strongest. Fig. 1131 shows that 
the role of antenna increases non-monotonically as a function of probability of 
spontaneous emotional arousal p s . We suppose that this phenomenon follows 
from a selection of positive emotions that the antenna receives from its neigh- 
bourhood. Such positive emotions are more likely for large value of p B where 
every agent fluctuates very fast. 

3.3 Stochastic resonance for emotionally driven task exe- 
cution efficiency 

It is generally accepted that emotions influence efficiency of various tasks 
executions, since the presence of negative emotions can sometimes make it im- 
possible to complete a difficult task while positive emotions can motivate a 
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Figure 15: Number NC of successful agents over a time T = 1024 with upper 
potential barrier U u = 30 and lower potential barrier Ud — — 50 as a function of 
p and p s for Model 2 and r = 2. 
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Agents efficiency in Model 1 
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Figure 16: Number NC of successful agents over a time T = 1024 with upper 
potential barrier U u = 30 and lower potential barrier Ud — —50 as a function of 
p for Model 1, for t = 2 and different values of p s . 
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Figure 17: Number NC of successful agents over a time T = 1024 with upper 
potential barrier U u = 30 and lower potential barrier Ud — —50 as a function of 
p for Model 2, for t = 2 and different values of p s . 
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Figure 18: Number NC of successful agents over a time T — 1024 with upper 
potential barrier U u = 30 and lower potential barrier Ud — —50 as a function of 
p s 6 [0.001; 1] for Model 1, for r = 2 and different values of p. 
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Agents efficiency in Model 2 



900 



800 



700 



600 



500 - 




400 



p=0 
p=0.001 
p=0.01 
p=0.05 

p=0.1 
p=0.15 

p=0.2 
p=0.25 

p=0.3 
p=0.35 

p=0.4 
p=0.45 

p=0.5 
p=0.55 

p=0.6 
p=0.65 

p=0.7 
p=0.75 

p=0.8 
p=0.85 

p=0.9 
p=0.95 
p=1 



Figure 19: Number NC of successful agents over a time T — 1024 with upper 
potential barrier U u = 30 and lower potential barrier Ud — —50 as a function of 
p s 6 [0.001; 1] for Model 2, for r = 2 and different values of p. 
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person to find a solution of a difficult problem. Of course, counter examples 
can also be presented, e.g. a high arousal connected with positive emotions can 
act as an obstacle for efficient work while a feeling of fear can magnify efforts 
to escape from a dangerous situation. 

Here, we examine a simple model of emotional influence on task execution 
efficiency as follows. Let the scalar variable Ui(t) describe the progress in task 
completed by an agent i at time t. Let us assume that this progress is only due 
to the presence of a temporary positive emotion while a temporary negative 
emotion leads to a partial damage of already completed work. To make our 
model as simple as possible let us consider the discrete dynamics of the variable 
Ui(t) in the form: 

u i (t) = u i (t-l)+e i (t) (4) 

As the initial condition of every agent we take Uj(0) = 0. To complete his 
task the agent needs to reach at certain time moment t the level Uj(t) > U u 
and such a successful event can be considered as a kind of jump over an upper 
potential barrier U u . On the other hand, when iii(t) < U4 (a lower potential 
barrier) the agent experiences a failure. A success or a failure means termination 
of agent's actions and causes a replacement of the agent i by a new agent j so 
during the simulation the number of agents is constant. The new agent starts 
his task with Uj (t) = 

To quantify the system efficiency we observe the number NC of successful 
agents that completed their tasks during some time window T. Results of 
corresponding simulations are presented at Fig. 14 and 15 for various values of 
parameters p and p s . One can observe that the dependences NC{p) (see at Fig. 
[TBI and IT71) and NC(p s ) (see at Fig. fTSland fTUI) are not always monotonic and 
for specific values of p and p s there are maximal numbers of successful agents. 
We can say that for very small values of p and p s agents are not effective at 
work. It can be interpreted as follows: a high level of emotions that occurs 
when the parameters p and p s are large disturbs people to work but when this 
level is too low individuals do not have enough motivation to undertake difficult 
tasks. The effect can be called an emotion-task resonance and it reminds the 
phenomenon of stochastic resonance when there is an optimal level of noise for 
signal transmission by a nonlinear system. 

4 Conclusions 

We have investigated two-dimensional models of affective agents. In the 
limit of small noise parameter p s describing the effect of spontaneous emotional 
arousal there are large amplitude collective oscillations of mean group emotion 
< e > (t) although the system dynamics does not contain any inertial part. As 
result the mean value of the group emotion is equal to zero. We suppose that 
the presence of oscillations and the lack of the spontaneous symmetry breaking 
are due to individual relaxation processes. The characteristic oscillations period 
increases with the value of the relaxation time r and the oscillations amplitude 
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increases with probability p of local affective interactions. For small values of 
relaxation times r we observed a stochastic resonance, i.e. the signal to noise 
ratio SNR is maximal for a non-zero p s parameter. The presence of emotional 
antenna can enhance positive or negative emotions and for the optimal transition 
probability p the antenna can change agents emotions at longer distances. The 
stochastic resonance has been also observed for the influence of emotions on task 
execution efficiency, i.e. there are optimal values of p and p s parameters when 
agents are most successful at completing their tasks. It suggests that neither 
completely non-emotional nor very emotional atmosphere should be arranged 
as working conditions. When the transfer of emotion between agents is too high 
(people are disturbed) or to low (they have no motivation) the task execution 
efficiency is very low. It is interesting that the optimal p s and p parameters 
values are similar for both resonance phenomena. 
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